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■ number states against multiphoton Jaynes-Cummings model (JCM). In particular, we inves- 

tigate atomic inversion, photon-number distribution, purity, quadrature squeezing, Mandel 
^ ^ Q parameter and Wigner function. We show that the quadrature squeezing for three-photon 

absorption case can exhibit revivals and collapses typical to those occurring in the atomic 
'<—'■.. 

^ ■ inversion for one-photon absorption case. Also we prove that for odd number absorption 

■ 

^ ' parameter there is a connection between the evolution of the atomic inversion and the cvo- 
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In this paper we discuss the quantum properties for superposition of squeezed displaced 



lution of the Wigner function at the origin in phase space. Furthermore, we show that the 
nonclassical states whose the Wigner functions values at the origins are negative will be al- 
ways nonclassical when they are evolving through the JCM with even absorption parameter. 
Also we demonstrate that various types of cat states can be generated via this system. 
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^ ■ I. INTRODUCTION 

^ . Jaynes-Cummings model (JCM) [l| is the simplest nontrivial example of two interacting quan- 

tum systems: a two-level atom and a single mode of the radiation field. In addition to its being 
exactly solvable, the physical system that it represents has become experimentally realizable with 
the Rydberg atom in high-Q microwave cavities, e.g. [2]. Such experimental realization has stimu- 



lated an intense research devoted to highlighting and genera 
this model has been extended tc) include, e.g., multiphoton [3 

and damping and dissipation [7]. Investigation of the JCM has produced many interesting effects, 



iziiig the original model. For example, 
, intensity-dependent coupling ^ 



which have no classica' 
quadrature squeezing 



analogue, such as revival-collapse phenomenon in the atomic inversion 
and sub-Poissonian statistics [lO|. Furthermore, JCM has been used as 



a source for the nonclassical states by means of the conditional measurement technique 



lit]. For 



more details related to JCM reader can consult On the other hand, evolution of states with 
JCM is an important topic, which can give insight into the occurrence of the nonclassical effects 
and how can be controlled based on the values of the interaction parameters, 
can mention, e.g., Shrodinger-cat states , [iSl], displaced number states 

coherent states 

Recently, we have introduced a general class of quantum states as a single mode vibration of 
electromagnetic field suddenly sq ueezed-plus-displaced by a collection of two displacements out 



n this respect we 
16l | and squeezed 



of phase with respect to each other 



2l|. These states (SSNDS) are defined as 



\a, r, n)e = X^[D{a) + eD{—a)]S{r)\n), 



(1) 



where is a normalization constant, D{a) and S{r) are displacement and squeeze operator, 
respectively, while a and r are displacement and squeeze parameters, e = |e| exp(i0), where (p 
is the relative phase and throughout the paper e = 0,i,l,— 1 in correspondence with squeezed 
displaced number state, Yurke-Stoler-, even- and odd-squeezed-displaced number states. Also |n) 
denotes Fock state. Squeeze and displacement operators are given, respectively, by 



S'(r) = exp[-(a^ — d^^)], Z)(a) = exp(a^a — aa*), (2) 
where d and d^ are annihilation and creation operators. The normalization constant is given by 



|A,|-i = l + |e|2 + 2|e|exp(-2|/?|2)L„(4|/3|2)cos</., (3) 

with 

P = a cosh r + a* sinh r ; (4) 

and are the Laguerre polynomials. Furthermore, the density matrix Pf{0) of this state can 

be written as 

pj{0) = \^f\0)){4^f\0)\ 

(5) 

= \Xe\'^{pM + Pi), 

where {0)) = \a,r,n)^ and the subscript / stands for the radiation field density matrix. The 
second line of ([5]) includes the statistical mixture part pM and interference part pj of the density 
matrix. The statistical mixture part has the form 

PM = D{a)S{r)\n){n\S^r)D^a) + \e\'^D{-a)S{r)\n){n\S^ {r)D^-a), (6) 
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while the quantum interference part is 



PI = e 



*D{a)S{r)\n){n\S\r)D\-a) + eD{-a)S{r)\n){n\S\r)D\a). 



(7) 



This part, i.e. pi, of the density matrix includes information about the quantum interference 
between the components of the state. On the other hand, in the number state basis ([T]) takes the 



form 



21 
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|a,r,n), = ^ Cm(a,r,n, e)|m), (8) 

m=0 

where the distribution coefficient has the form 

t h 

^' ^) = v21Z ^^ + (-1)^"+'"^^] exp[^^(tanhr - 1)] 

(9) 

min(m,n) ■\ f \ / \ 

y n!m! r 2 i y [ — tanhr i ^"^''■^ xj / 1 TT / e''a \ 

^ ^f^Q i!(n-i)!(m-i)! lv^er2;:J i 2 J ""'^-J' \^ ^^M^^ ^'"-^^ \^ ^/ilS^^ ' 

whereas Hm(-) is the Hermite polynomial of order m and a has been considered to be real. The 
photon-number distribution P{m) of the state ^ is 

P{m) = \Cm{a,r,n,e)\\ (10) 

For completeness, it has been shown that states ([T]) can be generated by the so-called quantum 

nn 

state engineering and also by trapped ion [20|, |21[ . The quantum properties of these states reveal 
that they can exhibit sub-Poissonian statistics, quadrature squeezing and oscillations in photon- 
number distribution. Their phase properties from the point of view of the Pegg-Barnett Hermitian 
pha^e fetalis, fl show ...»he..e»e two regimes controlling the behaviour based on whether 
the superposition is macroscopic (a >> 1) or microscopic {a < 1). For such superposition the 
nonclassical effects are more pronounced in the microscopic regime. Moreover, the influence of 



thermal noise on the behaviour of ([T]) has been considered in [23| showing that the correlation 
between different oscillators is essentially responsible for the occurrence of the nonclassical effects, 
which is similar to that of Schrodinger-cat states with thermal noise j24l ]. 

For convenience we mention that a new set of squeezed states using group-theoretical methods 



has been presented earlier, e.g. 25|. Actually, the definition of this set is based on the Holstein- 
Primakoff realization of both SU{2) and SU{1, 1). It has been shown that such type of states 
can exhibit interesting squeezing properties, depending in a characteristic way on the dimension of 
the irreducible unitary representation adopted. Moreover, the multiphoton, many-mode squeezed 
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states for SU{n), using a generalization Holstein-PrimakofF realization have been also reported [26l |. 
Finally, a long list of references related to different types of squeezed states can be found in the 
recent review \2^. 

In this paper we study the evolution of the states ([T]) with the multiphoton Jaynes-Cummings 
model. The motivation behind this work is to determine the evolution of the nonlinear process 
of the multiphoton transition in the presence of quantum interference of two squeezed displaced 
number states. Also, as we will show below, there are some intriguing features such as information 
about atomic inversion of the JCM can be obtained from quadrature squeezing of the nonlinear 
JCM, and also there is a connection between the evolution of the Wigner function at the origin 
in phase space and the atomic inversion. In addition to that this work gives a generalization to 
several results given in the literatures earlier, e.g. [l^-jiol. Actually, there is some merit in having 
the most general results. 

The Hamiltonian controlling the system in the rotating wave approximation (RWA) is 

^ =uJoa'a + uJa^z + gia^^+ + a^^cj^), (11) 

where a± and are the Pauli spin operators; loq and uja are the frequencies of cavity mode and 
the atomic transitions, respectively; g is the atom-field coupling constant and k is the absorption 
parameter. 

Throughout the paper we assume that the atom is initially prepared in its excited state | +) and 
the field is in the SSDNS ([T|). In this case the initial density matrix for the system is 

m = \+){+\6^\4\o))i^f\o)i (12) 



where stands for the direct product. Working in the interaction picture and considering the 
resonance case ujq = koJa-, after well-known analytical procedures, the time evolution density matrix 
describing the system is 

/)(r) = |^(r))(g)(v^(r)| 

oo f 

= E Cm{a, r, n, e)C^, (a, r, n, e) | cos[T^/h{m, k)]\m, +) (13) 

m,m'=0 

—i sin[T /i(m, k)]\m + k, — )| (^|(m', -|-| cos[T^/h{m', k)] + i{m' + A;, — | sm[T^/h{m', A;)]|, 

where T = gt is the scaled time, |— ) is atomic ground state and h{m, k) = -^-^r^- Also \'ip{T)) is the 
dynamical state vector of the system. From p3p the entanglement of the two interacting quantum 
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systems is readily apparent. To evaluate quantities associated with atom (or field) we have to trace 
(I13p over the field (or the atom) system. As it is obvious the system includes various parameters, 
which make such analysis difficult. For this reason throughout the paper we restrict ourselves to 
the parameters, which can give significant results. Also throughout the text the statement standard 
JCM means k = 1, the optical cavity field and atom are initially prepared in coherent states and 
in the atomic excited state, respectively. 

The paper is organized as follows: In section 2 atomic inversion is investigated. Section 3 is 
devoted to the photon-number distribution and purity. In section 4 Mandel Q parameter and 
quadrature squeezing are discussed. In section 5 Wigner function is demonstrated. Conclusions 
are summarized in section 6. 



II. ATOMIC INVERSION 



In this section we discuss the revival-collapse phenomenon in the evolution of the atomic in- 
version for the density matrix (jlSp . This phenomenon is a pure quantum mechanical effect and 
laving its origin in the granular structure of the photon-number distribution of the initial field 



281 1. Furthermore, this phenomenon has been realized experimentally in the sense that the state 
of the atomic beam leaving the cavity is monitored by ionization detectors [2]. On the other 
hand, the revival-collapse phenomenon has been seen also in nonlinear optics for the single-mode 



29( 1 when the modes are initially prepared in 



mean-photon number of the Kerr nonlinear coupler 
coherent light, however, in this case the origin of the occurrence of such phenomenon is in the 
presence of nonlinearity in the system (third-order nonlinearity specified by the cubic susceptibil- 
ity). We should stress here that the analytical formulae of both the mean-photon number of the 
Kerr coupler and of the atomic inversion of the JCM (under the harmonic approximation) [^] are 
similar. This remark makes the occurrence of the collapse-revival pattern in the Kerr coupler is an 
expected result. 



Now the atomic inversion for the density matrix (jlSp reads 

(a,(r)) = Tr/[a,(0)p(r)] 



(14) 



P{m) cos[2T y/h{m,k)], 

m=0 



where P[m) is given by (jlOp and the subscript / on the right hand side means that we trace 
over the field system. Firstly, it is worth reminding that the behaviour of the {(Jz{T)) of squeezed 
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FIG. 1: The atomic inversion (a) and the corresponding initial photon- number distribution (b) for the 
Yurke-Stoler displaced number states and for a = 5,/j = 1; the different values of n are specified in the 
figures. In Fig. la E denotes that this curve is associated with the even displaced number states whose 
P{m) is given in Fig. lb as dashed curve. Also ((T^(T)) and P{m) are shifted from bottom to top by 0, 2, 4, 6 
and 0, 0.09, 0.23, respectively. 

coherent light exhibits various interesting effects [itI . [l^. For instance, it has been shown that 
for strong coherent contribution [17] the cohapse time depends on the direction of the squeezing 
and for certain values of the squeezing parameter the response of the atom is similar to that of 
the chaotic radiation field. On the contrary, for strong squeezing contribution [3] the dynamical 
response of the atomic inversion shows echoes after each revival resulting from the interference 
effect. We start our discussion for k = 1 and the cavity field is in superposition of displaced number 
states (see Figs. 1 for given values of the parameters). Actually, the superposition of displaced 
number states can exhibit strong sub-Poissonian statistics as well as quadrature squeezing. Also 



they can be generated via the driven Jaynes-Cummings model or via trapped ion [301 . \3l\ . We 
proceed, it has been shown that there is a close correspondence between the number of peaks (with 
Poissonian envelopes) involving in P{m) and the number of occurred revivals in {az{T)) [16]. In 
other words, the envelope of each revival is a readout of the photon distribution, in particular, for 



slowly varying photon distribution 



321 ] provided that the atom is in the excited or in the ground 



state. Actually, for appropriate values of the parameters and when e = i, the connection between 
the peaks in P{m) and the revivals in {(Jz{T)) can be explained as follows [3]: The different 
peaks in P{m) have slightly different local mean photon numbers fhj (see Fig. lb) and thus each 
peak of P{m) gives its own revival-collapse pattern, which interfere with that of the others to 
produce a very complicated structure. In this case we have a sequence of the revivals with slightly 
different periods Tj = •n^frnj^ j = 0, 1, .... The echoes in the first revival will be separated by 
the intervals (Tj+i — Tj) whereas those in the second revival will be separated by 2(Tj_|_i — Tj) 
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FIG. 2: The atomic inversion {az{T)) for the asymptotic (curve 1) and the exact 2 + ((T^(T)) (curve 2) 
cases for even displaced number state with a — 16, fc = 1, r = and n = 1. 



and so on. Such a systematic behaviour depends upon the chosen values of the parameters, in 
particular, n (compare different curves in Fig. la). Furthermore, from this figure one can observe 
that the ringing revivals are remarkable for the case n = 2. For large values of n we noted that the 
chaotic behaviour is dominant and the collapse time is shortened. Now we turn our attention to 
the superimposed case, i.e. the even or odd case (see curve E in Fig. la). From this curve it is clear 
that the collapse time is smaller than that of the statistical mixture case owing to the interference 
in phase space. Precise information about this situation can be obtained from the asymptotic form 
for the atomic inversion, which can be evaluated via the harmonic approximation technique [181 ]. 
For this purpose and for the sake of simplicity we confine ourselves to|a| >>l,fc = l,n = l and 
r = 0. In this case the argument of cos(.) in ()14|) can be expressed as \XS\ : 



\/n + 1 = ■\/(n) + n + 1 - (n) ~ ^(?yi + rj2n), 



(15) 



where 



1 



(16) 



and (n) is the initial mean-photon number for the state under consideration. On substituting ()15p 
into (jl4p and after some minor algebra we arrive at 
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FIG. 3: The atomic inversion for even displaced Fock state (E-curve-((T2(r))) and odd displaced Fock state 
(O-curve-2 + (<t^(T))) for a = 5, fc = 2, n = 1 and r = 0. 



{a,{T)) = Ae{(l + |e|2)/i(r) cos[r(7?i + m) + |apsin(r/2T)] 
+2|e|/2(r) cos</>cos[r(r?i +772) - |ap sin(?72r)]|, 



(17) 



where 



/i(T)= [l-4|a|2sin2(2f) 



/2(r)= [l-4|apcos2(2^)' 



exp 



exp 



-2|a|2sin2(2^) 



-2|a|2cos2(2^) 



(18) 



In (jl7p . fi(T) and f2{T) are the envelope functions associated with the statistical-mixture part (jG]) 
and the interference part ([7]) of the optical cavity field density matrix, respectively. These functions 



are periodic with period 2'K^J (h) and their maxima occur for T^"^^ = 2t: yjln) and T)^' = TTyJin) 
where the superscripts {S) and (/) mean that the times are associated with the statistical-mixture 
and interference parts, respectively, and the subscript R denotes the revival time. It is obvious that 

1. Further, the term 1 



Td^^ = 7)T^\ which is typical with that of the Schrodinger-cat states case 



R ~ 2' 

in the square bracket of (jlSp corresponds to the coherent component in the optical cavity, whereas 
the second term is coming from the Fock state In other words, the main revival associated 
with the initial coherent light is splitted into three parts showing subsidiary revivals (echoes) in 



the atomic inversion, e.g. for the statistical-mixture part splitting occurs at sin( 



-2\a\- 
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FIG. 4: The atomic inversion (a) and the corresponding initial photon-number distribution (b) for a = 
3,e = i, k — I, n — I and different values of r as indicated. {(Tz{T)) and P{m) are shifted from bottom 
to top by 0, 2 and 0, 0.26, correspondingly. 

Furthermore, we have checked the case n = 2 and found that, e.g., the statistical-mixture-part 
envelope splits into four parts controlled by the equation 



I |4 ■ q| |2 ■ 2/%^^ , 1 

\a\ sm ( ) — 8|a| sm ( ) -|- 1 



0. 



(19) 



2 ' ' ' '2 

It is evident that the Fock state |2) develops two additional terms in the asymptotic form (jlOp . 
which provide more echoes in {az(T)) compared to that of the case |1). In Fig. 2 we compare the 
behaviour of the asymptotic formula (I17p with that of the exact one (I14p for the given values of the 
parameters. This figure reveals agreement and disagreement between the exact and asymptotic 
behaviour. Qualitatively, they are approximately in agreement in the locations of the revival- 
collapse pattern as well as each revival includes three portions, but the corresponding revivals 
cannot provide typical amplitude and shape. For the exact case (curve 2) the revivals are broader 
than those for the asyniptotic case. This problem can be overcome by means of the higher-order 
approximation in ()15|] [18t |. however, the price that should be paid, is that we cannot obtain a 
simple closed form for {az(T)). On the other hand, excluding the initial revival in Fig. 2, the 
origins of the first and second revivals, respectively, are statistical-mixture and interference parts 
of the initial density matrix of the radiation field. In conclusion the interference in phase space 
decreases the collapse time and the excitation number n yields echoes. 

The issue we would like to discuss here is that, qualitatively, the nature of quantum interference 
in phase space cannot manifest itself in the behaviour of the atomic inversion for the one-photon 
JCM. This is not the case for the multiphoton JCM where, e.g., the difference between the even- 
state and odd-state cases becomes clear, in particular, for k = 2. Actually, the two-photon JCM 
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has been extensively studied and extended in many direction (e.g. [33]). Information about the 
case = 2 is shown in Fig. 3 for given values of the interaction parameters. Clearly, the revivals 
in the two-photon case are much more compact than those in the one-photon case. Physically 
speaking in the two-photon process the phase correlation of the two superimposed optical cavities 



is transferred to the atomic dynamics 15|]. This is in contrast with the one-photon case where 
the atom absorbs only one photon at a time. Furthermore, such behaviour can be realized via the 

n 

asymptotic form of this case, where we use the asymptotic expansion |15l | 

3 



V(n + l)(n + 2) ^^+n. (20) 

Comparison between (jlSp and (j20p leads to the fact that the asymptotic results related to the 
two-photon case can be obtained from those of the one-photon case by simply setting 771 = 3 and 
rj2 = 2. In this case the envelope functions revive periodically with period T = vr and then the 
revival time is T = tt/2. So that the collapse and revival times become approximately intensity 
independent and this is responsible for such a systematic behaviour in the evolution of the atomic 
inversion. On the other hand, we have noted that for k > 3 the revivals interfere to produce chaotic 
population after the first collapse Q], i.e. the collapse-revival phenomenon of the atomic inversion 
for multiphoton Jaynes-Cummings model is more dramatic. 

We conclude this section by demonstrating that including squeezing in the superimposed optical 
cavity can stimulate revival-collapse pattern in the behaviour of {(Jz{T)). In doing so, we plot Figs. 
4 for the shown values of the parameters. From Fig. 4b one can observe that for r = 0, P{m) 
exhibits two peaks, which are smooth Gaussian peaks. Such structure leads to chaotic behaviour 
in {(Jz{T)) (see Fig. 4a for the corresponding case) and can be explained as follows: the second 
peak in P{m) gives the familiar collapse-revival pattern in {az{T)), whereas the first peak (which 
is incomplete one) provides oscillatory behaviour. Therefore, the competition between these two 
processes produces such behaviour. On the other hand, when r ^ 0, i.e. squeezing in the optical 
cavity starts to play a role, one can observe that P{m) has oscillatory behaviour with smooth 



Gaussian envelope, which causes the well-known revival 



situation is similar to that of squeezed number states 



16fi 



■collapse phenomenon in (o"z(T)). This 



III. PHOTON-NUMBER DISTRIBUTION AND PURITY 



In this section we discuss the evolution of the photon-number distribution P(m, T) and purity Tf 
for the system under consideration. As is well known photon-number distribution can be measured 
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FIG. 5: The photon-number distribution for e — 0,a — 7,k — 1 and n — 1 and for (a) T = 22.21; (b) 
T = 20.64. 

by photon detectors based on the photoelectric effect. Also, the characters of the field can be 
emphasized via P{m,T), which only at particular values of the interaction time resemble those 
of the initial pure states. On the other hand, investigation of the purity for the present quantum 
dynamical system is important to know how far the atom and the field are entangled. 
We start the discussion by analyzing P(m,T), which has the form 

P(m, T) = P{m) cos [Tv^/i(m, k)] + P{m - k) sm'^[Ty/h{m - k, k)], (21) 

where P[m) are given by (jlOl) . As the optical cavity field starts to interact with the atom the 
initial field distribution changes because of the quantum dynamics described by the Hamiltonian 
()lip . however, at specific times the field and the atom subsystems return to their initial states. 
Basically such behaviour depends upon the values of both the interaction time and the absorption 
parameter. To be more specific, this behaviour is established well for k = 2 and 4, as we will show 
shortly. On the other hand, for the standard JCM at one-half of the revival time pure-superposition 
states can be g enerated and the optical cavity field and the atom become completely disentangled 



36 



37( 1 . These states are well established in the strong-intensity regime, i.e. large values of 



|a|. Now we discuss the possible occurrence of such behaviour in the present system. To achieve 
our goal in a convenient way we restrict ourselves to the case e = and r = 0. In figures 5 we plot 
P(m, T) for T = vri/(n(0)) = 22.21 (one-half of the revival time) (a) and 7rA/(n(0)) - O.Svr = 20.64 
(close to one-half of the revival time) (b) and for the given values of the parameters. Actually, at 
these times the optical cavity field will be approximately in pure states this will be noticeable in 
the behaviour of the purity (see Fig. 7 below). Furthermore, Figs. 5 reveal that the optical cavity 
collapses to superposition of displaced number states, particularly, superposition of D{a)\l). This 
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is obvious where P{m 
to such type of states 



T) has oscillatory behaviour with two-Gaussian envelopes representative 



30( 1 ■ Of course the oscillations in P{m,T) are resulting from the quantum 



interference in phase space. The quantum properties for such type of states (,i-e., superposition 



of displaced number states) have been discussed in [30[ and they have been generated also in a 



state-selective measurement [3l|. We proceed, comparison of Fig. 5a and Fig. 5b reveals that the 
oscillations in P{m, T) are extremely sensitive to the values of the interaction time. Furthermore, 
comparison between the oscillatory behaviour in the two peaks in Fig. 5a (or in Fig. 5b) shows 
that a clear asymmetry and this is connecting with the asymmetry in the initial photon-number 
distribution of the states D{a)\l) (see Fig. lb). Actually, by changing the types of the initial state 
we can obtain different forms of cat states. 

Now we discuss the case k = 2 (two-photon absorption) whose scaled revival time is Tr = vr, 
where we still consider e = 0. We found that the optical cavity field at T = sTr where s is 
positive integer, collapses to pure state typical (or similar to) the initial states, while the atom 
becomes either in the ground state or in the excited state. In other words, the atom and the 
field subsystems become completely disentangled at these specified times. Such behaviour can be 
analytically determined provided that a is large and n is finite through the substitution of both 
T = Tr = vr and ()20p into the state vector obtained from (jl3p . we arrive at: 



\^{T = Tr))= Z Cm(a,r,n) cos[ri?(m+|)]|m,+) -«sin[rjj(m+f)]|m + 2,-) 

m=0 ^ 

(22) 



= |Vf (0))(g)exp(ivr)|-), 

(2) 

where (0)) denotes the field state having the form 

oo 

l^fiO)) = Yl C„^{a,r,n)exp{^m7^)\m + 2). (23) 

m=0 

Expression ()23p shows that the cavity field possesses photon distribution as that of the initial one, 
but with two-photon shift. On the other hand, by applying the above procedures for the case 
T = 2Tr = 2tt, one can easily show that the system returns back to its initial form, i.e. the atom 
is in its excited state and the field is in \ip^p (0)) given by ([1]) with e = 0. Generally, we conclude 
that the system evolves periodically according to the following relation: 



(U I IV'I^^(O)) for T = sTr and s = 1,3,5,.., 

\^?in = { ' (24) 
I IV'/ (0)) iorT = sTR and s = 2,4,6,..,. 
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FIG. 6: The photon-number distribution for k — 2, T — jTt and the other parameters as those in Fig. 5. 



Expression (j24p can be understood as follows: when the interaction is going on the atom decays 
(spontaneously and stimulately) from the excited state to the ground state contributing its photon 
to the cavity field and the system becomes disentangled during interaction time T = Tr. This 
means that the atomic excitation energy is transferred to the field 36|], i.e. an information about 



the initial state of the atom is stored in the field. Furthermore, when the interaction continues the 
feed back between field and atom occurs and after T = 2Tr the system returns back to its initial 
form. This behaviour is also visible in the evolution of the atomic inversion and reflects the periodic 
nature of the system. We should stress that the origin of such behaviour is two-fold: (i) We have 
considered the system is completely isolated (the interaction with the environment is neglected), 
(ii) We treat the system in the RWA, which ensures that whenever k photons are lost from the 
field the atomic state must change from ground to excited states, or vice versa. Furthermore, there 
is a similarity between the present behaviour and that of the frequency conversion device in which 
the energy interchanges periodically between the signal and idler modes. On the other hand, it 
has been shown that for the case k = 2 and at the quarter of the revival time superposition states 
can be generated [sfi]. Presently we show analytically using a similar technique as that in 34] that 



similar behaviour can occur. We restrict ourselves to the case e = 0, n = 1 and a is very large. On 
substituting T = jTr = jir and (j20p into the state vector (|13p . we obtain 
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n = 



n=l ° 




k=3 



k=2 



FIG. 7: The purity against the scaled time T for: a) (a,fc, e) = (7,1,0) ; and b) (a,n, e) = (7,1,0). In (a) 
for n = we plot Tf + 0.5. In (b) for fc = 3 and 4 we plot Tf + 0.5 and Tf + 1, respectively. 



\^P{T = \Ti^)) = Y: C^{a,r,n)\cos[\TR{m + ^)]\m,+) -ism[\TR{ 



m+ 4)l|m + 2 



(25) 



cos[jTR{m + |)]|Cm(a,r,n)|+) + iCm-2{a,r,n)\-)\ 



m) 



m=0 



For the case n = 1 the relation between Cm{o!,r,n) and Cm-2{o!,r,n) reads 

2' 



Cm-2{a,r,n) 



m 



m — a 



m — 



C.m{a,r, n). 



(26) 



a a 

It is worth remembering that the initial photon distribution has two-peak structure (see Fig. lb) 
each of which has a Poissonian envelope. Further, for very large values of a one can find out that 
the terms contribute effectively to the summation in (j25p are those for which the values of and 
m are comparable. Therefore, the first term in the square root of (|26|) tends to unity while the 
second one tends to zero. Also the terms in the square bracket reduces to unity where » 2. 
Considering all these suggestions, ([25]) can be rewritten as 



1 OO T "i 

mr = -Tr)) = J2Cmia,r,n) cos[-TR{m + -)]\m) (g){|+) + (27) 

m 

From this expression it is clear that at the quarter of the revival time the atom and the field become 
asymptotically disentangled and the atom evolves in an atomic superposition state, whereas the 
cavity field collapses to superposition of states whose forms are somewhat similar to those given in 
33]. From ([27l) it is clear that cos[jTR{m + |)] 7^ for m integer meaning that the P{m,T) has 
oscillatory behaviour but cannot have perfect oscillations (,i.e., P{m,T = Tr/A) = 0) more than 
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those included in the initial distribution of the optical cavity. This behaviour is clear in Fig. 6 where 
we plot the P{m,T) at T = Tr,/4 for the same situation as in Fig. 5 but k = 2. Furthermore, in 
Fig. 6 the solid and dashed curves are associated with the exact ([21]) and asymptotic (f27ll forms of 
P{m, T). From Fig. 6 one can observe that the solid and dashed curves include significantly similar 
behaviour. In other words, under these conditions the asymptotic expression (I27p can be used to 
describe the optical cavity field well. It is worth mentioning that approximate perfect oscillations 
can be achieved in the behaviour of this case if one uses the replacement Tr Tr it 0.08, however, 
the behaviour of P{m, T) will be very close to that in Fig. 5b. Furthermore, in (jlip if replaced 
by aVa^ (i.e., the system becomes intensity-dependent JCM) the system can generate perfectly 
even and odd superposition of displaced number states (considering the initial state is the displaced 
Fock state) according to the atom is prepared initially in the ground and in the excited states, 
respectively. 

We close this section by analyzing the purity for the optical cavity field via the relation Tf = 
Tr[/5j(T)] where Pf{T) = Tra[/5(T)] is the field reduced density matrix, whereas p{T) is the total 
density matrix of the system (I13|) . When Tf = \ the reduced density matrix describes pure state, 
however, for Tf < 1 the field will be in a statistical mixture state. The smaller the purity is, 
the less pure the state will be. Further, when Tf = 0.5 the field and the atom will be in a 
maximally entangled state. As we mentioned above Tf gives information about entanglement of the 
components of the system, which forms the basis of experiments in the realm quantum information. 
Further, the disentanglement of the two quantum systems suggests interesting applications, e.g., in 
atomic states preparation through interacting quantum systems or in the choice of optimum time 
of flight in a micromaser to determine the cavity field 39(]. Moreover, from the point of view of 
the quantum theory of measurement, the decoupled field state can be regarded as a pointeer basis 



since all the information about the atomic state has been transferred to the field [3j]. We proceed, 
in Figs. 7 we plot the evolution of the Tf for the given values of the interaction parameters. 
Comparison between the curves associated with n = 1 and n = in Fig. 7a (for fixed value of 
k) shows that including Fock state |1) in the coherent optical cavity give rises to the oscillatory 



behaviour in Tf, changes the parity and achieve maximal entanglement at T = 2^^^ (^(0)) (revival 
time). Actually, the oscillatory behaviour in Tf occurs in a good correspondence with the revival 
in {az{T)). Further, for both cases (n = 0, 1) the field collapses to a pure state at the middle of 



the collapse time {t^^/ (n(0))). Also the figure does not show a well-defined ringing revival for the 
case n = 1, which is in contrast with the behaviour of the {(Jz{T)). Furthermore, we have noted 
that increasing the value of n will not give significant change in the overall behaviour of Tf except 
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increasing the oscillations. However, for the superimposed optical cavity, we have noted also that 
the state of the field never gets as pure as in the Fig. 7a and the system remains in a maximally 
(or partially) entangled state most of the time. On the other hand, Fig. 7b is given for fixed value 
of n and different values of k as indicated. For the case k = 2, Tf evolves smoothly and the field 
collapses to an approximate pure state at T = ir/4 and 37r/4, which is long lived as compared 
to those occurring at times multiples of vr. As we mentioned above for T = 311,3 = 0, 1,2.., the 
optical cavity field returns to its initial form. Furthermore, the field and the atom are long lived 
maximally entangled at T = tt/2 and 37r/2. It is evident that Tf has a systematic behaviour as 
{az{T)). Now we turn our attention to the case k = 4. One can observe that Tf displays pure states 
instantaneously with period it/2, however, at T = 7r/4 and 3it/4, the field will be in a partially 
mixed state. Apart from these values of interaction time the field and the atom are approximately 
maximally entangled. The similarity and dissimilarity between the Tf for the cases k = 2 and k = 4 
are remarkable. On the other hand, for the case k = 3, the field goes to approximate mixture states 
just after the interaction is turned on meaning that the field and the atom almost retain a strong 
entanglement all over interaction times. 

In conclusion, the behaviour of P{m, T) and Tf are in a good agreement with each others. Also, 
the analysis of the purity gives information that could not be predicted by looking at the inversion or 
the photon-number distribution solely. We have not considered the infiuence of squeezing parameter 
in the present section since it generally increases the oscillatory behaviour in the P{m,T). 



40( 1 having the 



IV. SUB-POISSONIAN STATISTICS AND QUADRATURE SQUEEZING 

In this section we discuss the sub-Poissonian statistics and quadrature squeezing for the system 
under consideration. Both of these quantities can be used as a measure for the nonclassical effects. 
Firstly, the sub-Poissonian statistics can be analyzed using Mandel Q parameter 
form 

^(^. _ ((An(T))^) - (n(T)) 

Formula (j28p characterizes the deviation from Poissonian statistics. It holds that — 1 < Q{T) < 
for sub-Poissonian statistics (nonclassical effect), Q{T) > for super- Poissonian statistics (classical 
effect) and Q{T) = for Poissonian statistics (the standard one). Sub-Poissonian light can be 
measured by means of a set of photodetectors. Further, this light has several applications, e.g. in 
the gravitational wave detector and quantum nondemolition measurement, and can be generated 
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FIG. 8: Mandel Q parameter against \a\ when the optical cavity field is initially prepared in D{a)\l) (a) and 
even coherent state (b) for T = 1.578 and for k = 1, 2, 3, 4 corresponding to solid, short-dashed, long-dashed 
and star-centered curves. In (a) the bell-centered curve is for zero-interaction time. Straight line is the 
sub-Poissonian bound. 



in semiconductor lasers 411] and in the microwave region using masers operating in the microscopic 
regime 

Secondly, in order to investigate squeezing we define two quadrature operators as X = 
^[a + a^, Y = ■^[a — a'^]. These quadratures relate to the conjugate electric and magnetic field op- 
erators of the electromagnetic wave and satisfy the commutation rule [X, Y] = ^. Therefore, the un- 
certainty relation reads {{AX{T))^){{AY{T)f) > ^. So we can say that the system is able to pro- 
duce squeezing if the squeezing factor F(r) = A{{AX{T))^)-1 < or S{T) = A{{AY{T))^)-1 < 0. 
As is well known squeezed light can be measured by a homodyne detection in which the signal is 
superimposed on a strong coherent beam of the local oscillator. Furthermore, quite recently it has 
)een shown experimentally that there is an evidence of squeezed light in the biological systems 
43|. 

The required expectation values to treat these quantities can be calculated using the following 
relation 



{a^'iT)a''{T)) = Tr [p(r)at«(0)a^'(0) 



(29) 



where p{T) is given by (jl3p and s,s' are positive integers. Before discussing the sub-Poissonian 
statistics for the present system it is worth reminding that for the standard JCM the sub-Poissonian 
statistics have been realized jlol |. Furthermore, the establishment of such statistics for input 
squeezed coherent states relies on the competition between the squeezing parameter r and co- 
herent amplitudes |a| [Si]. In other words, as the values of r increase the amount of sub-Poissonian 
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statistics (resulting from the coherent component) is decreased and eventuahy disappeared. For 
the system under consideration the squeezing parameter plays the same role, i.e. it destroys the 
sub-Poissonian statistics. In general, we found that Q{T) exhibits collapse-revival pattern as well 
as ringing revivals, which reflects the time evolution of the atomic inversion. In this case the photon 
statistics already oscillates between sub-Poissonian and super-Poissonian statistics based on the 
values of the interaction parameters. Here we would like to address two facts: (i) The system can 
enhance the sub-Poissonian interval for the initial sub-Poissonian states, (ii) The system can gen- 
erate sub-Poissonian light for an initial non-sub-Poissonian superimposed state (e.g. even coherent 
states). Information about these facts is included in Figs. 8a and b for given values of the parame- 
ters. In these figures we have plotted Mandel Q parameter in the stationary regime at specific value 
of interaction time, since the time-dependent coefficients in the present system are periodic in time. 
The selected value of the interaction time gives maximum sub-Poissonian statistics for the case 
k = 1 when the optical cavity field is initially prepared in the even coherent states. From Fig. 8a, 
one can observe that the curves are smooth and the most sub-Poissonian statistics occur for k = 3 
and 4. Also this figure gives the coherent amplitude interval over which the sub-Poissonian light 
can be obtained for different values of the absorption parameter. It is evident that the intervals 
for the cases k = 1,2,4 are greater than those for the initial one (compare different curves with 
the bell-centered one). This shows that how one can enhance the sub-Poissonian interval based 
on the value of the absorption parameter. Furthermore, the amount of sub-Poissonian statistics 
decreases almost exponentially as the magnitude \a\ increases similar to that of the standard JCM 



4J]. Figure (b) carries information about the behaviour of the Mandel Q parameter when the 
optical cavity field is initially in even coherent states. It is worth reminding that this state exhibits 
always super-Poissonian statistics. However, here through the evolution we can observe that the 
sub-Poissonian statistics are achieved giving their maximum values for A; = 1,4. In Fig. 8b the 
curves include smooth oscillations resulting from the superposition principle. Further, the sub- 
Poissonian intervals are increased compared to those in the Fig. 8a, in particular, for k = 1,2,4. 
The two-photon absorption case (short-dashed curve) almost exhibits steady-state sub-Poissonian 
statistics for large values of \a\, whereas the case A; = 3 exhibits sub-Poissonian statistics only for 
the large values of |a|. Finally, figures 8 show that there is no clear relation between the obtained 
sub-Poissonian statistics and the values of the absorption parameter k. 

Now we turn our attention to discuss quadrature squeezing. Actually, it has been shown that 
squeezed light can be obtained in the standard JCM Q]. However, for input squeezed light with the 
atom initially prepared in excited atomic state the interaction destroys the initial squeezing, which 
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(b) 




FIG. 9: The squeezing factor F against \a\ (a) and scaled time T (b) for even coherent states. In (a) 
{k,T) = (1,12.65) (solid curve), (2,4.5) (short-dashed curve) (3,7.7) (long-dashed curve), (4,9.6) (star- 
centered curve) and the straight line is the squeezing bound. In (b) (fc, a) — (3, 6) and the inset is the 
atomic inversion of the case (fc, a) — (1, 6). 



is inherent in the radiation field [19|. Here we note that squeezing can be obtained provided that 
the initial superimposed optical cavity is squeezed. Moreover, the system can switch squeezing from 
one quadrature to the other. The origin of this is in the periodic nature of the system. For example, 



even coherent state with real a exhibits squeezing in Y-quadrature provided that < |q| < 2 
In the present system squeezing can be transferred to X-quadrature independent of the values of 
k (see Fig. 9a for given values of the parameters). In this figure the values of the interaction time 
have been chosen numerically in order to maximize the amount of squeezing in F. We proceed, 
from Fig. 9a for each value of k, F starts from almost minimum-uncertainty state around |a| = 
to increase the amount of squeezing until it exhibits its maximum squeezing, then starts to decrease 
the squeezing and eventually exhibits classical effect (i.e. squeezing completely disappears). Also 
from this figure the maximum (minimum) amount of squeezing occurs for the case k = 2 (1) 
for these chosen values of the parameters. Fig. 9b includes the dynamical behaviour for the X- 
quadrature for the case (A;, \a\) = (3,6). It is obvious that there is no nonclassical negative values 
(, i.e. there is no squeezing). Nevertheless, F{T) exhibits collapses and revivals, which (apart 
from the different scales) are very similar to those of the ((Tz(T)) for the case k = 1 (compare 
the behaviour of the inset in Fig. 9b with that of F{T)). To be more specific, qualitatively both 
these curves are equivalent in such a way that they collapse, remain quiescent, revive, collapse 
again and repeatedly undergo a complicated pattern of collapses and revivals for large interaction 
time, which is not included in the figure. Such type of equivalence leads to an interesting result: 
the collapse-revival phenomenon in the {az{T)) for the standard JCM can be detected indirectly 
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via measuring the quadrature squeezing for three-photon absorption case through the homodyne 
detection. It is worth mentioning that observation of collapses and revivals were performed using 
the one-atom mazer j^, which is more sophisticated than the dynamics of the JCM. 

V. WIGNER FUNCTION 



In this section we investigate the behaviour of the Wigner (1^) function for the system under 
consideration. Actually, the W function is informative, sensitive to the interference in phase space, 
and can give a prediction to the possible occurrence of the nonclassical effects. Furthermore, this 
function can be constructed experimentally using the optical homodyne tomography j46|. 



The W function is defined 



47[ as 



W{x,p,T) = — / dCexpHpC)(x+^|/5/(r)|x-|), 



(30) 



where we have assumed that h = ujq = 1. On substituting the wavefunction obtained from (|13p 
into (jSOp and using both the wavefunction of the Fock state as 



{x\n) 



V7r22"n! 



(31) 



and the identity 

/oo 
exp(-C')H^(C + zi)H„(C + Z2)dC = 2^V^m\z^-"'L^-"'{-2z,Z2), m < 
-oo 



n, 



(32) 



we arrive at 



W{x,p,T) = '-^^stmi g (-l)™'x"'2^C™(a,r,n,e)C;;,(a,r,n,e) 

m,m'=0 



cos 



iTy^hi^] cos[ryR^]yglL--™'(2|x| 



(33) 



where x = 2; + ^p. 

Before discussing the behaviour of (j33p it is worth mentioning that the Q distribution function 
of the standard JCM with and without cavity damping have been investigated in 49|, ISOl . ISll . l52l |. 



It has been shown, for the non-damping case, that the initial shifted Gaussian distribution splits 
into two distributions, which counter-rotate on a circle in the complex plane of the distribution. At 
the opposite end of the circle the two peaks collide, they split again, and so forth. This behaviour 
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is closely connected to the revival-collapse phenomenon in the {(Jz{T)), where collapses of the 
Rabi oscillations occur during splitting of the distribution and the revivals occur during collision 
of the peaks of the distribution. As the system evolves the width of the peaks become broader. 
When they spread over the whole circle no collision occurs anymore and distinct collapses and 
revivals are no longer observed. Also splitting in the Q distribution for the JCM is explained 
based on the interpretations of the density matrix 53]- Further, it is worth mentioning that the 



generation of superposition states in the standard JCM has been investigated using this function. 
For the damping case the counter-rotating peaks spiral into the origin till finally the stationary 
distribution is reached. Employing W function for JCM has stimulated few efforts, in particular. 



for the damping case 



49[. Here we provide a connection between the evolution of both {az{T)) 



and the value of the W function at the origin for the two following cases. 

Case (i): when the absorption parameter is odd number, i.e. k = 1,3,5, ... From (133]) we can 
write the following expression 



Vr(0,r) = - y (-Ij'^PM cos[2T ^Jh{m,k)], (34) 

m=0 

where only the diagonal terms survive. Careful examination to (j34|) leads to the following facts: 
for the pair-photon states (i.e. even or odd states) ([M|) can be rewritten as 



M^±(o,r) = ±-(<T,(r)), (35) 

vr 

where the positive (negative) sign refers to even (odd) case. This expression indicates that 
when W{Q,T) ~ (7^ 0) the collapse (revival) occurs in {(Jz{T)). It is obvious that when 
VF(0, T) = zbl/vr the amplitude of {az{T)) is maximum. Furthermore, the value of W function 
at the origin exhibits nonclassical negative values in the course of revival times, which means 
that the nonclassical effects most probable occur during this time. On the other hand, for 
the non-pair-photon states the opposite behaviour occurs. For instance, for the standard JCM 
VF(0, T) 7^ (—0) the collapse (revival) occurs in {az{T)). These ensure that Schrodinger-cat 
state can be generated during the collapse time, where the establishment of the interference in 
phase space is quite obvious. This fact will be made clear shortly. 

Case (ii): when the absorption parameter is even number, i.e. k = 2,4,.... In this case the 
origin of the W function is localized in phase space (time independent) having its initial value 
W{0, T) = W{0, 0). This leads to a novel result, which is: the nonclassical states whose W functions 
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FIG. 10: W function of the even displaced Fock state ior a — 3, k = l,n = 1 and for: a) T = 4.599998 
(collapse time) and h) T — 12.60001 (revival time). 



have negative values at the origins (e.g. odd coherent state) will be always nonclassical when they 
are evolving through the JCM with even absorption parameter. Of course the chaotic behaviour 
(and/or overlap of the successive revivals) in {(Jz{T)) can be recognized in the evolution of the W 
function at the origin when |Ty(0, Ti)| ~ |Ty(0, T2)| such that Ti > T2. It is worth referring to the 
nonlinear version scheme of a single-atom homodyne detector in which the ionization probability 
of a single test atom is given in terms of the Wigner characteristic function of the field [5j]. In 
Figs. 10 a and b we plot W function for even displaced number state for T = 4.599998 (collapse 
time) and 12.60001 (revival time), respectively, for the given values of the parameter. The values 
of the interaction time are chosen numerically. Actually, W function of one-photon even coherent 
state consists of two inverted holes, with more pronounced negative values, and interference fringes 
inbetween. In Figs. 10 according to the evolution of the dynamical system such shape is destroyed 
and the W function becomes more structuralized, which is more pronounced for the collapsed case 
Fig. 10a. To be more specific. Fig. 10 a includes multi-peak structure ensuring the generation of 
the Schrodinger-cat states. Also for this case we find that W{0, T) ~ 0. Conversely, Fig. 10 b 
includes more pronounced negative values at the origin of phase space. These facts are in a good 
agreement with the above conclusions. 

The second issue we would like to address here is that it has been shown that cat states can 
be generated at one-half of the revival time in the standard JCM 



341, 



35 



36 



33]. One of the 



ways is the investigation of the Q function. Nevertheless, this function does not include a complete 
information about the system as the W function includes. Now we give the form of the W function 
for such type of cat states. We do so for the standard JCM, i.e. k = 1 and e = n = r = 0, with 
strong intensity in such a way that the harmonic approximation is applicable. To deal with the 



23 



second part of ([33|) some tricks should be performed, e.g. replacing Cm by Cm+i since we are 
working in the strong- intensity regime such a transformation is relevant. On using (llSp in addition 
to the properties of the Laguerre polynomial (f33ll reduces to 



W{x,p,T) = ^exp[-(3;-ei)2]|exp[-(p-e2)^] +exp[-(p + e2)^] 

(36) 

-2exp(-p2)sin(^)sin [T(r/i - f ) + 2a (y2x - acos(^)) sin(^)] }, 



where 



ei = V2acos(^), £2 = \/2asin(^). (37) 

We should stress that for the best of our knoweldge this is the first time for which the exact form 
of the W function, i.e. (|36p . of the generated cat state in the microcavity is given. Expression (j36p 
shows that the peaks centers move on the circle controlled by (|37p . In fact, the appearance of the 
factor -v/2 in the definition of the circle is related to the definition of W function (j30p . however, if 
one uses the transformation {x,p) —>■ V^{q,y) then the relations (p6]l and (f37|l would be consistent 
with those obtained earlier for the Q function, e.g. [52], provided that the prefactor l/(27r) is 
replaced by l/vr. Such a prefactor is connected with the transformation in the integration (j30p . 
Expression (j36p shows that at the revival times (i.e. T = 2-7r/?72) the two peaks collide, whereas at 
one-half revival time the well-known shape of the cat-state W function is obtained, which is two 
Gaussian bells corresponding to statistical mixture of individual composite states and interference 
fringes inbetween arising from the contribution of the overlapping between different components 
of the state. It is worth referring to |52], in which an analytical closed form for the Q function 
is derived by applying the method of steepest decent and saddle-point integration in the limit of 
high photon numbers. The final remark, the asymptotic form for the W function related to the 
two-photon absorption case can be obtained from (j36p by simply setting rji = 3 and 7^2 = 2. In this 
case the W function is intensity independent and its behaviour is systematic as that of {az(T)). 



VI. CONCLUSIONS 



In this paper we have discussed the quantum properties for superposition of squeezed displaced 
number states defined by ([T|) against the lossless multiphoton Jaynes-Cummings model. We have 
assumed that the atom is initially prepared in the excited state. In addition to the characteristic 
features of the system the investigation gives generalization to various results in the literatures. 
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The investigated quantities are the atomic inversion, photon-number distribution, purity, Mandel 
Q parameter, quadrature squeezing and W function. The main conclusions are the followings: 
Atomic inversion exhibits revival-collapse pattern as well as echoes, which are readout of the 
photon-number distribution. Furthermore, including squeezing in the superimposed optical cavity 
can cause revivals and collapses in the atomic inversion. The photon-number distribution shows 
that various forms of cat states can be generated via the system. Also the entanglement and 
disentaglement between the field and atom occur periodically since the system is characterized by 
the periodic exchange of energy between the atom and the cavity mode. Such behaviour has been 
confirmed in the purity evolution. Actually, the behaviour of these quantities basically depends 
upon the values of the absorption parameter k, which results in the functional dependence of the 
Rabi frequency on the mean photon number. The sub-Poissonian statistics occur regardless of the 
values of k provided that a is finite. Also squeezing in the quadrature squeezing can be occurred 
and switched between the two quadratures provided that the initial superimposed optical cavity 
is squeezed. Furthermore, the evolution of the quadrature squeezing for the case A; = 3 exhibits 
revivals and collapses similar to those of the atomic inversion of the standard JCM. Actually, this is 
a novel result and its consequence is that the revival-collapse phenomenon exhibited for the atomic 
inversion can be measured by homodyne detection via three-photon microcavity [3]. With respect 
to W function we have found some novel results, which can be summarized as follows: There is a 
direct relation between the evolution of the values of the W function at the origin in phase space 
and the evolution of {cFz{T)) when the absorption parameter k is odd number. This suggests that 



the atomic inversion can be detected via both the homodyne tomography 



and trapped ion 



5a]. It is worthwhile mentioning that in the trapped ion methods a proposal for measuring the 
W function at the origin of the phase space for a single photon field is given [56;] and this would 
be relevant for the task. On the other hand, for even absorption parameter the origin of the W 
function in phase space is localized and thus the nonclassical states whose W functions include 
negative values at the origin will be always nonclassical when they are evolving with the JCM. 
Furthermore, for the first time, as far as we know, we have given the explicit form of the W function 
of the generated cat state in the microcavity. 
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